clear
set more off

use "E:\Cuong - paper\COVID\Pvalue_governance_paper.dta"

	gen n=_n

	qqvalue q1_pvalue, method(simes) qvalue(q1_qvalue)
	qqvalue q2_pvalue, method(simes) qvalue(q2_qvalue)
	
	* qqvalue pvalue, method(yekutieli) qvalue(qvalue1)

	br

	graph  hbar (asis) q1_pvalue q1_qvalue, over(name) bargap(5) legend( label(1 "P-value") label(2 "Q-value") ) ///
	ytitle("Value") title("P-value and Q-value of the lowest quintile") ylabel(0(0.1)1) blabel(bar, format(%4.2f))
	
	
	graph  hbar (asis) q2_pvalue q2_qvalue, over(name) bargap(5) legend( label(1 "P-value") label(2 "Q-value") ) ///
	ytitle("Value") title("P-value and Q-value of the near lowest quintile") ylabel(0(0.1)1) blabel(bar, format(%4.2f))
	
	